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Abstract. The response of a neo-Hookean fiber composite undergoing finite out-of-plane sliear 
deformation is examined. To tliis end an explicit close form solution for the out-of-plane shear 
response of a cylindrical composite element is introduced. We find that the overall response of 
the cylindrical composite element can be characterized by a fictitious homogeneous neo-Hookean 
material. Accordingly, this macroscopic response is identical to the response of a composite cylinder 
assemblage. The expression for the effective shear modulus of the composite cylinder assemblage 
is identical to the corresponding expression in the limit of small deformation elasticity, and hence 
also to the expression for the Hashin-Shtrikman bounds on the out-of-plane shear modulus. 

keywords: fiber composites, hyperelastic composites, nonlinear composites, finite elas- 
ticity, effective properties, micromeclianics 



1. INTRODUCTION 

The problem of characterizing the behavior of fiber reinforced solids undergoing large deforma- 
tions arises in many modern applications such as biomechanics and active materials. A desirable 
goal is to express the behavior of these materials in terms of the properties and the volume frac- 
tions of the constituents composing the composite. This will allow, for example, to account for 
variations in the volume fractions of the phases and to obtain estimates for the stresses developing 
in the constituents. Motivated by this goal, in this work we study the problem of characterizing 
the behavior of fiber reinforced composites under out-of-plane shear loading. This type of loading 
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may arise, for instance, during a torsion of a cylinder made out of a composite reinforced with 
circumferential fibers. 

Following the work of Hashin and Rosen (1964), who considered the analogous problem in the 
limit of small deformation elasticity, we examine the problem of a single cylindrical composite 
element {e.g., Hill, 1964) subjected to shear loading in the fiber direction. For the case of a 
composite element made out of two incompressible neo-Hookean phases we provide a simple close 
form solution describing the strains and the stresses in the two phases. These expressions correctly 
reduce to the small deformation limit and are in agreement with the results of Hashin and Rosen 
(1964). The solution can be easily extended to coated composite elements, that is the case where 
the fibers are coated with a few concentric shells made out of different materials. 

We further demonstrate that, under this type of loading, the behavior of the composite ele- 
ment is identical to that of a fictitious homogeneous neo-Hookean material. An expression for the 
effective shear modulus of this fictitious material, in terms of the properties and the volume frac- 
tions of the phases composing the composite element, is determined. This expression is identical 
to the corresponding expression of Hashin and Rosen (1964) for the out-of-plane shear modulus 
of a composite cylinder assemblage in the small deformation limit. By following arguments simi- 
lar to those of Hashin and Rosen (1964), we also conclude that this shear modulus describes the 
macroscopic behavior of a composite cylinder assemblage undergoing finite out-of-plane shear de- 
formations. Finally, we recall that in the limit of small deformation elasticity the expression for the 
effective out-of-plane shear modulus is identical to the expression for the Hashin-Shtrikman bounds 
of Hashin (1965) on the out-of-plane shear modulus of two phase transversely isotropic composites. 

2. OUT-OF-PLANE SHEAR DEFORMATION OF A CYLINDRICAL COMPOSITE ELEMENT 

We consider a cylindrical composite element made out of a fiber with a circular cross section 
surrounded by a concentric shell of a different material. The elastic behaviors of both phases are 
characterized by neo-Hookean strain energy-density functions 

W^'\A) = i/iWTV(A^A - I), (1) 

where s = 1 for the fiber phase and s = 2 for the shell. The radius of the fiber-shell interface is a, 
the outer radius of the shell is b, and we assume a perfect bonding at the interface between the two 
phases. It is assumed that the diameter of the cylindrical element is at least an order of magnitude 
smaller than its length. We choose a coordinate system such that in the reference configuration the 
X3 axis is aligned with the direction of the fiber. The composite element is subjected to a simple 
out-of-plane shear deformation such that in the deformed configuration the location of the points 
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on the outer boundary of the sheh is given by the mapping 

Xi = Xi, X2 = X2, X3= X3+ 7X1, 



(2) 



where 7 is the amount of shear (Ogden, 1997). In a cyhndrical coordinate system, where the Z 
axis is ahgned with the fiber direction and the ray = with the positive direction of the Xi axis, 
the above mapping is expressed in the form 



i?, 61 = 6, z = Z + jRcosQ, 



(3) 



where {R,Q,Z) and {r,9,z) are the coordinates of the boundary points in the reference and the 
deformed configurations, respectively. 

For the material points in the two phases we consider the following mappings 



(4) 



where the functions f^^^ are functions of i?^*^ only. The continuity condition across the interface 
resulting from the third of Eq. (4) is 

/«(«) = /(2)(a), (5) 

and at the outer boundary of the shell the boundary condition results in the requirement 

/(2)(6) = 6. (6) 

In cylindrical coordinates the expressions for the deformation gradients in the two phases are 
{e.g., Humphrey, 2002), 

^1 o\ 

A(^) =0 10 (7) 

^7/{;'^ cos e(^) -i7/(*)sine(*) ly 

where //j is the derivative of / with respect to its argument. We note that det A^'^^ = 1. 

In terms of A'-'^^ the nominal (or Piola-Kirchhoff) stresses in the two neo-Hookean phases are 

/ ^i^)-pi^) 7pW/R^coseW \ 

S(") = /xW-p(^) -i7pW/WsineW (8) 

V7/^^'^/^^cose(^) -i7;uW/WsinGW ^iW-p^^^) y 

where n^^^ is the shear modulus of the s-phase and p^^^ is the arbitrary pressure in the phase. In 
accordance with the assumption that the length of the cylindrical composite element is substantially 
longer than its diameter, far from the bases the pressure is independent of the Z coordinate, that 
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is p*-*^ = p^^\R, 0). The traction continuity conditions across the interface between the two phases 



are 



^(i)-p«(a,e) = /i(2)-p(2)(«,e), 

(g) 

and we note that the traction continuity condition for Srq is satisfied. 

For convenience, we follow Ogden (1997) and express the equilibrium equations in terms of the 
nominal stress. Assuming no body forces, in cylindrical coordinates {e.g., Humphrey, 2002) these 
equations lead to the following three equations, namely, 

— — = — = ' H 1 — = (10) 

OR ' R de ' i?2/ ^ R dR ^ dR^ ^ ' 

From the equilibrium along the R and @ directions it follows that the pressure in each phase is 

constant, and from the equilibrium along the Z direction 



/W(i?) = ci^)i? + r^. (11) 

We complete the solution by noting that in the fiber phase we must choose = 0. The 
remaining three integration constants in Eq. (11) are determined from Eqs. (5), (6) and the second 
of (9). Accordingly, in the fiber phase 

2;u(2) 

Z = Z + ^- 7— 7— -— 7-7- i? COS 6, (12) 

'(1 -c{i))^(i) + (l + c(i))^{2) ^ ^ 

and in the shell 

7^ (MW+^(^))-cW(^)^(/^W-/i(^)) . ^ 

^ = ^ + ^ (i-ca))^a) + (i + caV2) ^''^ 

where c^^^ = (a/b)'^ is the volume fraction of the fiber phase in the composite element. The jump in 
the pressure is determined from the first of Eq. (9). To determine the absolute values of the pressure 
additional information concerning the stresses must be specified. Finally, the explicit expressions 
for the stresses in the two phases can be readily determined. 

3. COMPARISON WITH A FINITE ELEMENT SIMULATION 

We compare the analytic solution developed in the previous section with corresponding numerical 
finite element simulation. Taking into account symmetry conditions, a quarter of the cylindrical 
composite element (0 < © < 7r/2) was constructed and meshed by application of the commercial 
finite element code ABAQUS. The ratio of the length of the cylindrical composite element to its 
diameter was 5. The volume fraction of the fiber was c^^^ = 0.25 corresponding to a/b = 2. The 
nodes in the {R, Z)-plane containing the ray = were constrained from moving in the Q direction. 
The nodes in the (i?, Z)-plane containing the ray Q = tt/2 were constrained from moving in the G 
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Figure 1. Out-of-plane displacements (a) and shear stresses (b) as functions of R 
along the rays = and 7r/4. The continuous curves correspond to the analytic 
solution and the marks to the finite element simulation. 



and the Z directions. The boundary condition according to Eq. (2) for the circumferential nodes 
were setup manually for 7 = 0.4. This value of 7 corresponds to a principal stretch ratio A = 1.22 
{e.g., Ogden, 1997). ABAQUS built in neo-Hookean constitutive model was used and the ratio 
between the shear moduli of the fiber phase and the shell was /x^^^ //^^^^ = 10. 

Shown in Fig. la are the out-of-plane displacements Uz = z — Z as functions of R along the rays 
G = and 7r/4. In Fig. lb the variations of the out-of-plane components of the Cauchy stress Trz 
and Tqz are shown as functions of R along the same rays. In both figures the continuous curves 
correspond to the analytic solution and the marks to the numerical results determined at the nodes 
of the finite element mesh. To minimize edge effects, the numerical results were extracted at the 
central cross section of the model half way between its bases. For convenience, the displacements 
are normalized by the out-of-plane displacement Uz{b, 0) of the circumferential point along the ray 
= 0. The stresses are normalized by the uniform stress component Trl\r,0) developing in the 
fiber along the ray = 0. 

Clearly, there is a fine agreement between the results from the analytic solution and those from 
the numerical simulation. In the numerical simulations results we observe the small jump in the 
stress component T0z{r,O) at the interface between the fiber and the shell. We verified that this 
jump is mesh sensitive and diminishes as the mesh becomes finer. 
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4. Applications to fiber reinforced composites 

An important question concerns the ability to identify an effective homogeneous medium whose 
behavior is identical to the behavior of the cylindrical composite element. To examine this issue 
we follow Hashin and Rosen (1964) and determine the traction developing at the circumferential 
boundary of the cylindrical composite element under the deformation (4). From Eq. (8) for the 
nominal stress the components of the traction on the boundary are 

tR = ^i^^^ - p^^\ te = 0, tz = J fi cose, (14) 

where 

,,, (l + cW)^W + (l-cW)^(^) 

^ ^ (l-cW)//(l) + (l+c(l))/i(2)- ^""^ 

Consider next a homogeneous cylinder made out of a neo-Hookean material with shear modulus 
^0- The cylinder is subjected to the same boundary conditions as given in Eq. (3). The corre- 
sponding deformation gradient and the stress in the cylinder are constants. It can be easily verified 
that the components of the traction developing on the boundary of this cylinder are 

tR = fJ'0-Po, te = 0, tz=^fiocosQ, (16) 

where po is an arbitrary pressure. 

We note that if = /^j up to a constant arbitrary pressure, the overall out-of-plane shear 
response of the composite cylinder and the homogeneous cylinder are identical. Accordingly, the 
out-of-plane shear behavior of a composite cylinder assemblage (Hashin and Rosen, 1964) made 
out of a space filling assemblage of cylindrical composite elements with various radii but with fixed 
ratio of the inner to the outer radii in each element will be identical to the shear response of a 
homogeneous neo-Hookean material with shear modulus fl. For the idealized transversely isotropic 
composite cylinder assemblage the above homogenized result for the out-of-plane shear modulus is 
exact. For more realistic microstructures with overall transverse isotropy Eq. (15) can be used as 
a straightforward and simple approximation for the out-of-plane shear modulus. 

5. Concluding remarks 

Within the framework of finite deformation elasticity, we develop a close form solution for the 
problem of a cylindrical composite element made out of two incompressible neo-Hookean phases 
subjected to out-of-plane shear loads. We note that since the shear problem is isochoric, if the 
behavior of one of the two phases in the cylindrical composite element is neo-Hookean and the 
behavior of the second phase reduces to a form identical to that of a neo-Hookean material under 
shear loads (even if, in general, the second phase is compressible), the proposed solution will satisfy 
this problem too. This is because under the imposed boundary conditions the deformation of the 
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second phase must be isochoric. We further note that the proposed solution is in agreement with 
the corresponding solution of Hashin and Rosen (1964) for the analogous problem in the limit of 
small deformation elasticity. We also demonstrate that the results of a corresponding finite element 
simulation are in agreement with the proposed solution. 

In the context of fiber composites undergoing large deformations, we find that the effective out- 
of-plane shear behavior of a composite cylinder assemblage (Hashin and Rosen, 1964) made out of 
two neo-Hookean phases is identical to that of a homogeneous neo-Hookean material. An expression 
for the effective shear modulus of the homogeneous neo-Hookean material in terms of the properties 
of the two phases and their volume fractions was determined. This expression is identical to the 
expression obtained by Hashin and Rosen (1964) for the analogous problem in the limit of small 
deformation elasticity. We also note that this expression is identical to the expression obtained by 
Hashin (1965) for the bounds on the effective out-of-plane shear modulus of transversely isotropic 
composites in the limit of small deformation elasticity. 

Finally, it is noted that this expression is also identical to an exact expression that was recently 
obtained by deBotton (2005) for the effective shear modulus of a special class of incompressible 
transversely isotropic composites undergoing finite in-plane shear deformations. In the limit of small 
deformation elasticity, the microstructure considered by deBotton (2005) is an optimal microstruc- 
ture attaining the Hashin-Shtrikman bounds of Hashin (1965) on the in-plane shear modulus for the 
class of incompressible transversely isotropic composites. Thus, we observe an interesting analogy 
where in two works dealing with transversely isotropic composites undergoing finite in-plane and 
out-of-plane shear deformations, the expressions for the effective shear moduli are identical to each 
other and to the ones previously obtained in the limit of small deformation elasticity. 
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